Two-Channel Kondo Physics from Tunnelhng Impurities with Triangular Symmetry. 



Aris L. Moustakas and Daniel S. Fisher 
Department of Physics, Harvard University, Cambridge, MA 02138 

Tunnelling impurities in metals have been known for some time to have the potential for exhibiting 
Kondo-like physics. However previous models based on an impurity hopping between two equivalent 
positions have run into trouble due to the existence of relevant operators that drive the system away 
from the non-Fermi-liquid Kondo fixed point. In the case of an impurity hopping among positions 
with higher symmetry, such as triangular symmetry, it is shown here that the non-Fermi-liquid 
behavior at low temperatures can be generic. Using various bosonization techniques, the fixed point 
is shown to be stable. However, unlike the conventional two-channel Kondo (2CK) model, it has four 
leading irrelevant operators, implying that while the form of the singular temperature dependence 
of physical quantities is similar to the 2CK model, there will not be simple universal amplitude 
ratios. The phase diagram of this system is analyzed and a critical manifold is found to separate the 
. non-Fermi-liquid from a conventional Fermi liquid fixed point. Generalization to higher symmetries, 

such as cubic, and the possibility of physical realizations with dynamic Jahn- Teller impurities is 
discussed. 
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. I. INTRODUCTION 

(N ■ 

Following pioneering work by Nozieres and Blandiro, the multichannel Kondo model has attracted much attention 
^ because of the prospect of novel, low temperature behavior. Various theoretical methods have been developed to study 
QQ : the natuffijof the zero temperature Kowlp-fixed points for a variety of cases, in particular with an arbitrary number of 
channelsa~Ll and with either one or twou Kondo impurities. Despite various intriguing predictions, no experimental 
PsJ , realization of the multichannel Kondo effect has yet been conclusively demonstrated. The difficulty lies in the fact that 
C ■ ' the non-Fermi-liquid fixed points are unstable to various symmetry-breaking processes, such as channel anisotropy, 
which turn out to be present in real experimental situations. For example, Nozieres and Blandiro pointed out that 
anisotropy between the channels, causedJay lattice effects, would destroy the non-Fermi-liquid ground state. 

Subsequently, Vladar and Zawadowskit3 were the first to suggest a different realization of the two-channel Kondo 
. model, with the role of the spins and the channels of the conventional magnetic model interchanged. In this case, a 
non-magnetic impurity tunnelling between two sites plays the role of the magnetic spin-1/2 impurity in the original 
^ I Kondo formulation with electrons coupled orbitally to the impurity, while the spin index of the electrons is a spectator, 
' playing the role of the channel index. Therefore, since in the absence of a magnetic field in the magnetic Kondo system 
^ ] the spin- up andj^in-down electrons are degenerate, channel anisotropy is no longer a problem. This led Vladar 
O ■ and Zawadowskillj to predict that such a system should exhibit non-Fermi-liquid behavior at low temperatures, 
' characteristic of the two-channel Kondo (2CK) model. Recently, low temperatup. tunnelling data by Buhrman et 
^ . a£3~tj in very small metallic contacts have been interpreted by Buhrman et atMB in terms of two-channel Kondo- 
like physics. Their measutements are claimed to be consistent with certain exac t Eesp lts for the 2CK model obtained 
by Affleck and LudwigJ3i3 but at this point, the interpretation is controversial Jl30 Again, the problem arises from 
' the absence of a symmetry which prevents other relevant terms from appearing; the most important of these is bare 
impurity tunnelling, which in the conventional Kondo model corresponds to a magnetic field on the impurity. In 
a previous paperEZI, we studied this system in detail and argued that it is very unlikely that 2CK physics could be 
observable over a substantial temperature range in such a two-site impurity system. 

A way to circumvent the problem of the appearance of relevant operators is to impose symmetry conditions on the 
orbital degrees of freedom so that the relevant operators which plagued the two-site impurity system are not allowed 
to appear by symmetry. The first question that naturally arises is what are the minimum symmetry requirements 
at some high (or intermediate) energy scale so that the system will fiow under renormalization (RG) towards the 
non- Fermi liquid fixed point, without any fine tuning. Clearly, the low temperature non- Fermi liquid behavior one 
should be primarily aiming for, is that of the two-channel Kondo model due to its relative simplicity. 

For tunnelling non-magnetic impurities, the "channel" (i.e. spin) symmetry is ensured, in the absence of a magnetic 
field, due to time reversal invariance. As for the symmetry of the impurity motion, an essential requirement is that 
at intermediate energies the impurity, dressed with high energy electrons, should behave like a degenerate doublet, 
just as the impurity spin-i does in the conventional 2CK model. For this to be true without any parameter tuning, 
the symmetry group of the system has to allow for degenerate doublets to occur, i.e. it has to have a two dimensional 
representation. The simplest group with such a property is the triangular group C31,, which, with six elements, is the 
smallest non-abelian group. 
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In this paper we are thus led to consider a non-magnetic impurity in a metal that can tunnel between the corners 
of an equilateral triangle. In this case, the left and right moving impurity states are related by symmetry and are 
thus degenerate. Not surprisingly, in order for the system to be in the basin of attraction of the 2CK point, we must 
require that at some intermediate energy, the left-right doublet of the impurity is the lowest energy orbital state of the 
impurity. Naively one might expect that this is impossible, since the ground state of a particle in an external potential 
is very generally non-degenerate. However, in a solid the impurity interacts with other electrons, including localized 
atomic electrons, whose state is affected by the impurity motion. By taking these interaction effects into account, we 
explicitly demonstrate in some simple and potentially physically relevant examples, that it is indeed possible for the 
ground state of the impurity-plus-high-energy-electrons complex to be degenerate. 

We then analyze the weak electron-impurity interaction limit, determine the criterion for flowing under the renor- 
malization group (RG) toward the 2CK fixed point, and find the energy scale below which the low lying impurity 
doublet becomes substantially separated from the excited symmetric orbital state. Below this energy scale the im- 
purity complex can be regarded as a degenerate doublet, and the system can be mapped to a 2CK model with only 
irrelevant perturbations away from the intermediate coupling 2CK fixed point. As a result, the 2CK fixed point is 
accessible at low temperatures for generic values of the bare couplings provided that the conditions for an impurity 
doublet obtain. 

In order to demonstrate the existence of novel 2CK-like behavior at low temperatures, one must do more than 
show that the RG flows at zero temperature go to the 2CK fixed point. As shown by Affleck, Ludwig and others, 
the singular low temperature behavior of various physical quantities in the 2CK model is controlled by the leading 
irrelevant operators about the 2CK fixed point: i.e. by the flow towards the fixed point. We must thus show that 
the structure of the irrelevant operators for the tunnelling impurity system is similar to that of the 2CK model. 
Indeed, we will show that the scaling dimension of the leading irrelevant operators, which characterize the form of the 
singularities in physical quantities at low temperatures, is the same in the two cases. However the number of leading 
irrelevant operators is larger for the tunnelling impurity system. Analyzing the behavior close to the 2CK fixed point, 
we find four leading irrelevant operators, in contrast to the pure 2CK model, which has only two, due to the fact that 
our problem does not have the full symmetry of the conventional 2CK problem, U{l)charge x SU{2)orbitai x SU{2)spin- 
We thus conclude that since different physical quantities can have singularities controlled by different combinations 
of the various leading irrelevant operators, various amplitude ratios and scaling functions will be less universal than 
in the 2CK problem, although the form of the singular low temperature behavior will be the same. 

In real systems, there will of course generally be some-hopefuUy small-symmetry-breaking perturbations. By 
studying the relevant symmetry breaking perturbations, in particular strains, we conclude that their effective strength 
will generally be screened by the fluctuations at the 2CK fixed point, making their effects less important than wojuM 
be guessed by naively comparing energy scales. However, for the experiments by Buhrman and collaborators£3t3 
simple order-of-magnitude estimates show that the strains have to be extremely small for the 2CK fixed point to be 
approached. 

Within the space of full triangular symmetry, we find that by tuning the splitting between the impurity doublet 
and the singlet states, (e.g. by pressure) there will be a phase transition at zero temperature from the 2CK behavior 
to conventional Fermi liquid behavior. This transition is conjectured to be controlled by a fixed point with SU (3) 
orbital symmetry. 

Impurity configurations with high symmetry, such as cubic, are also discussed briefly, and found to be characterized, 
up to additional leading irrelevant operators, by the 2CK fixed point. 



A. Outline 



In the next Section we motivate and introduce the model we will be using. In Section [II , we discuss the importance 



of the impurity degeneracy and analyze three physically relevant situations, for which the impurity ground state can 
indeed be two- fold degenerate. The weak coupling analysis is performed in Section [V, at which the crossover from 



the 3-site problem to a degenerate 2-level impurity system, as well as Csu-symmetry breaking processes are discussed. 
In Section 0, the intermediate coupling regime for the 2CK system is studied, while in Section VI conjectures about 
the behavior on and near the the critical manifold separating 2CK from conventional Fermi liquid behavior are made. 
Finally, we conclude in Section VH with the potential relevance, o£, the model to experimental situations and, in 
particular, to recent experiments by Buhrman and collaborators.E3~c£l 
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II. MODEL 



In this Section we introduce an effective Hamiltonian for the impurity and conduction electron system, vahd at 
some intermediate energy scale substantially less than the conduction electron bandwidth. The effects of higher 
energy electronic processes have already been taken into account by introducing renormalized parameters in the 
effective Hamiltonian. It is crucial that at this intermediate energy scale all processes that are allowed by symmetry 
will be present, even if they were absent in the "bare" Hamiltonian. Therefore we have to include all processes that 
are relevant (in the RG sense) at either weak or intermediate coupling, since their amplitudes, albeit small initially, 
can grow under the RG flows. [Other processes that are irrelevant for both weak and intermediate coupling will end 
up yielding at worst non- universal order one corrections to the various quantities of interest.] 

Since the impurity-electron interactions are predominantly local, the electronic degrees of freedom which interact 
most strongly with the impurity at low energies can be viewed as centered around the three minimum energy positions 
of the impurity. Therefore for low enough energies we only need to consider three electronic wavefunctions for 
each energy and each spin. These can be represented as a symmetric state c^^.^ with spin a —H and energy e, 
transforming according to the trivial representation Fi of 6*3^, and an orbital doublet {c]^„^, c^o-ei' 
dimensional representation, T3. Neglecting weaker interactions of the impurity with other electronic Fs-doublets, Fi- 
singlets and antisymmetric F2-singlets will not change the behavior of the system. In an analogous way the impurity 
states can be described by a Fa-doublet d^, and the symmetric state d^. 

Before introducing the effective Hamiltonian 7i, let us summarize its symmetries. Impurity and electron numbers 
are separately conserved, d e^'^d and c e*^c. In addition, in the absence of spin-orbit coupling or magnetic 
fields, our system has spin-S'[/(2) and time reversal symmetry. The latter acts also on the orbital degrees of freedom: 
i — > —i; d^ <-*■ da; cl <-> cr. Finally, Ti. has the C31, symmetry, which includes exchange symmetry around an axis 
(cl ^ cn; di, ^ dii) and rotational symmetry by ±27r/3: cl exp [±27ri/3] c^; cu exp [^27ri/3] cj^; cs ^ cs 
and similarly for the d's. 

The Hamiltonian Ti. includes several types of terms: 

Ti. = Hel + Hint + Tihop + Timix (1) 

with Tiei the electron kinetic energy in the absence of the impurity 

/de 
c[„,CLae + c'j^a^CRa, + cl^^CSae , (2) 

with the energy e measured from the Fermi surface. Note that for convenience we consider a linear conduction band 
with a cutoff of order the bandwidth W. Although we will generally work in the S,L,R representation, in order to 
see the physical origin of various terms, it is convenient to first work in a basis of wave-functions centered around the 
three equivalent impurity sites, with 

c„ = (c5 + e-'("-i)^CL + e'("-i)^Cfl) (3) 

for n = 1, 2, 3. If the impurity is at rest, then the only allowed one electron terms are 

Vi ^ dldnci^Cna (4) 

and 

V2 T.'^ndnclc,,. (5) 

The other such terms consistent with the symmetries can be rewritten, using dl^dn = 1, in terms of these two and 
terms which mix the electronic states centered at the different sites: 

which must, in any case, be present due to the non-orthogonality of the ci^2,3 statesJll [We ignore the unimportant 
term ctc„]. 
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The bare impurity hopping term can be written 

^hI - f E 4'^. = I (24^5 - d^n. dldj) (7) 

where the di,2,3 and dL.R.s are related as in Eq(|). The sign of £ is of crucial importance. In particular, if all electron- 
impurity interaction terms are much smaller than 5, then the sign of £ determines whether the impurity ground state 
is an orbital singlet {£ < 0) or a Fs doublet (£ > 0). At energy scales such that the renormalized electronic bandwidth 
is smaller than \£\, in which case tunnelling events of the impurity from the ground state configuration to excited 
states will be suppressed, the low temperature behavior clearly depends crucially on the symmetry of the impurity 
ground state in the absence of interactions. 

Of course, as was noted in the Introduction, naively the sign £ should generally be negative to ensure a non- 
degenerate ground state. However, situations can exist for which the interaction of the impurity with atomic (high 
energy) electrons conspires to make the effective low energy impurity hopping matrix element £ positive and hence 
favor the doublet impurity ground state. Some illustrative examples are discussed in the next Section. 

In addition to the bare impurity hopping, various "electron assisted" impurity hopping terms (Tij^op) can occur. 

These may involve either terms of the form dl^dnCjCj or d^^dnclcj with i 7^ j. At the Fermi surface there are six 
independent such terms with coefficients Ai...Aq. In the L,R,S basis, these give rise to several terms that are 
diagonal in this basis: 

Vdiag = TtQi ^ (c'L^CLa - Cr^CRo^ {d\dL ~ d\dR^ (8) 

(T 

+ E (^cLcSa - C^RaCB^a " c[^cl„) {^d^ds - d\dR - d\di}j 

a 

+ T^Qs (cl^csa + cia^^La + Cr^crct^ (sd^dg - d^'j^dR - d\dj}j 
where cj^g., for fj, — L, R, S, are the "local" electron creation operators, represented as cjj^^ = c^^^ {x = 0) where 

clA^)^ j ^e'"'4<.. (9) 



effect of the first two of these terms is to 
Their amplitudes Qi, Q2, Q3 correspond 



are the incoming (outgoing) electronic s-wave states for x > {x < 0). T 
screen transitions between the different states L, R and S of the impurity, 
to linear combinations of various electron-assisted tunnelling processes {Ai . . . Aq) and, Qi and Q2 will be assumed 
to be positive. If Qi < 0, then it will be seen later that the system behaves like a ferromagnetic Kondo model [see 
discussion after Eq(^0|)]. In order for these important terms to be of appreciable magnitude (in Section ^ we will 
allow Qi to be of order unity), we focus on impurities with relatively fast (renormalized) tunnelling rates. 

The interaction terms Vi and V2 give rise to a number of processes in the L, R, S basis which change the state of 
the impurity; to be consistent with C^^ symmetry the total orbital "angular momentum" - with L = — l,i? = +l and 
S = - oi each of these terms must be zero (mod 3). The electron assisted hopping terms {Ai . . . Aq) also contribute 
to Voff-diag- We have 

Voff-diag = 27rAi ^ {d\dRC^j^^CLcr + h.c}j (10) 

(T 

+ 27rA2 ^ {d\dsc\^CLa + (i^^dsC^^CRcr + h.c. 
a 

+ 27rA3 ^ {d\dR(^^^csa + d\dRc}g^CR„ h.c. 

a 

-I- 27rA4 ^ {d\dsC^R„CSa + ^^^C^^CSo- + /i.c.^ 

a 

+ 27rA5 ^ [d^hdsC^Lcr^Rcj + d^RdsC^Rcr^-Lcj + /i-c.) 
a 

In the absence of impurity motion, Vdiag is zero, but Voj f-diag has contributions from both electron-impurity inter- 
actions and electron-assisted hopping and is thus still non-zero. 
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Of all the terms we have so far, A3 and A5 are the oajy processes, that break the 0(2)-symnietry in the ra-manifold 
(which, using different notation, was discussed in RefO in the context of the two-site tunnelling impurity). However, 
these terms are clearly Csv invariant and their effects will be discussed in later Sections. 

The terms we have included so far can be conveniently grouped into 

and the Vdiag and Voff-diag combinations exhibited above. The electron mixing term, as it is diagonal in the i, R, S 
basis, will not play a substantial role (in contrast with the two-site case of RefO). In addition to the various impurity 
one-electron coupling terms, the impurity can also affect electron-electron interactions, and two-electron assisted 

hopping can take place. Although these terms are strongly irrelevant for weak coupling as discussed in Ref.EZI, the 

(2) 

latter can play an important role near the intermediate coupling 2CK fixed point; we denote them as n^^^. Putting 
all the important terms together, we primarily study the effective Hamiltonian 

H =Ho + Vdiag + Voff-diag + T^'hop- (12) 

At this point, one can immediately guess which terms will play dominant roles at very low energy scales where, 
with £ positive, the higher energy symmetric impurity state dg will be strongly suppressed. In this limit, the most 
important terms will be Qi which is analogous to the Jz coupling in the Kondo system, and Ai which is analogous 
to the J± coupling. The only other one-electron term which does not involve ds is A3 which involves the cs electrons 
and breaks the 0(2) symmetry of the L,R subspace. 

In general there can also be Cs^-symmetry breaking processes, but, by assumption, these will be small. Nonetheless 
we will need to investigate their role. The most important terms correspond to bare impurity transitions from state 
L to R, such as 

(^MdldR + M*d^^dL) , (13) 

or from L and R to S, such as c?^ [d^ + dji) + + dj^^ ^5. [Note that d^d^ — dj^rfj? and (cZl — dji) + {d\ — d\^)ds 

cannot appear if the system is time-reversal invariant]. These will originate from various effects, including strains, 
other nearby static impurities etc. Of the above processes, the one corresponding to M will be the most important 
to our discussion, because it will persist even on energy scales below which £ becomes larger than the renormalized 
bandwidth and transitions to the S'-state are effectively quenched. The real part of M will be a linear combination of 
two independent bare impurity terms that are even under parity <-^- d^, i.e. an energy bias term 2d\di — d\d2 — d\d^ 
and an asymmetry in the hopping amplitudes d^d^ -\-d\d2., while in the imaginary part of M similar terms will appear 
which are odd under parity. Thus, two independent parameters have to be tuned for \M\ to vanish. Finally, note that 
asymmetries in the electronic states due to other impurities etc., can generate M under renormalization. 

III. JAHN-TELLER EFFECT, BERRY'S PHASE AND EXCITED IMPURITY CONFIGURATIONS: 

EFFECTS ON IMPURITY DEGENERACY 

As already mentioned, the degeneracy of the lowest lying effective impurity state at intermediate energy scales is 
crucial to our analysis. The generic case of a non-degenerate impurity ground state, corresponding to £ < 0, leads to 
uninteresting impurity screening and hence to ordinary Fermi-liquid behavior at low temperatures. Instead, a doubly 
degenerate r3 orbital impurity state will be our focus, and in this Section we will show how the necessary condition, 
£ > 0, is possible. First, we will analyze in two model cases the effects on the impurity motion of high energy atomic 
or valence electrons in solids, ignoring coupling to the conduction electrons. The impurity will be assumed to move 
slowly on the scale of the atomic electronic time scales. We can thus consider the adiabatic evolution of the electronic 
ground state as the impurity moves. In the second part of this Section we will demonstrate how, when we include the 
conduction electrons, the existence of intermediate excited impurity configurations between the low lying states d\^ 
d\, d\, can serve as a mechanism for changing the sign of the effective £ at lower energies. 

A. Effects of atomic electrons on impurity degeneracy 

We focus first on an orbital- and spin-singlet interstitial impurity, tunnelling between three equivalent positions-1, 
2 and 3-on the corners of an equilateral triangle. The underlying physical reason for the desired effect on the sign of 
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£ is the appearance of degeneracies in the states of the surrounding medium (electrons and neighboring nuclei) when 
the impurity is located at a position of special symmetry-in our case, the center of the triangle. Even if the impurity 
is forbidden energetically from occupying this central symmetric position, e.g. by the presence of a host lattice atom, 
the effect is manifested in the behavior of the wavefunctions of the system when the impurity is away from the center. 

The simplest example is an atomic electronic state on the impurity, which, with the impurity located in the center 
of the triangle, is doubly degenerate, i.e. an doublet. When the impurity is at, say, site 1, the degeneracy is 
broken. However, by adiabatically moving the impurity once around the triangle and back to site 1, the electron 
returns to its original state, but with-its wavefunction multiplied by a factor of —1. The extra phase acquired by the 
electronic state is just Berry's phase;u2l it will change the sign of the effective inipurity-plus-electron tunnelling matrix 
elements. 

With one electron in an appropriate energy level on the impurity, this effect can occur but the impurity will be 
magnetic. What is needed instead is an orbital (non-magnetic) electronic doublet on the impurity when it is at the 
center of the triangle. |— ■ 

This is a situation-which is a classical example of the Jahn-Tellem effect and was also proposed recently in this 
context by Gogolin.u Due to the coupling of the electronic doublet of the impurity and the neighboring atoms, the 
triangle-center position of the impurity is unstable to displacements of the impurity from its equilibrium position: this 
linearly splits the electronic degeneracy, but only costs elastic energy that is quadratic in the distortion. A massive 
impurity will be displaced to one of three minimum energy equilibrium positions on the corners of an equilateral 
triangle. But, since the impurity wavefunctions at these three positions will generally have a finite overlap, the 
impurity can tunnel between them. The phases of the tunnelling matrix elements add up to tt on tunnelling around 
the triangle, yielding an extra minus sign from the Berry's phase. This ensures that the tunnelling Fs-doublet will 
be the ground state of the impurity complex. This restoration of the degeneracy via tunnelling of the impurity is 
a dynamic Jahn- Teller effect-the triangular symmetry is not broken in the ground state and the original electronic 
doublet has merely been modified into an impurity orbital doublet. 

Unfortunately, orbital electronic doublet ground states at high symmetry positions of an impurity are not too 
common. Instead, in many cases the minimum energy configuration has magnetic character, due to Hund's rules 
splittings. Nevertheless, for large enough local lattice distortions due to the impurity, the effective tunnelling impurity 
ground state can still be a doublet. To explicitly demonstrate this effect, we consider a simple situation with four 
atomic electrons on the impurity, perturbed by the surrounding medium, with spin-orbit interactions neglected. For 
simplicity we consider only three electronic levels which, with the impurity in the central point with triangular 
symmetry, will have the symmetry of the L, R and S states, corresponding to angular momentum of m — ±1, (mod 
3) respectively. We denote the creation operators for these ]^ g- If the environment is such as to make the atomic 
one-electron S state the lowest by energy e^, two electrons will fill this and the other two will go in the degenerate L, 
R states. The simplest situation is with effective electron-electron interactions smaller than the splitting Ea- In this 
case, Hund's rules will give rise to ferromagnetic exchange, J, between the two L ox R electrons, resulting in a spin 
triplet ground state, with orbital angular momentum (mod 3), m=0, a doublet pair of spin singlet, m = ±1 
states (I^T-^i > ^"^^ l-^T-^i ^ higher energy symmetric singlet, ^S, excited state. Couphng to the one-electron 

S states will mix, e.g., I^ii?! — SiRi^ > in with the \L^Li > state with some amplitude K. These processes can be 
represented by a Hamiltonian for the impurity at the triangle center, 

Ti-centef^ -^a ^ CL^Sa'^Sa + J {a\^aRT^a^j^^aLi + h.c}j (14) 

+K ^atj^OLf alj^aLx + a^g-^aRi;a\^aRi + a^j^^aLia^g^ULi + a^i^^aRi;a^g^an.i + h.c 

[Note that if all four electrons are considered together, with the effective interactions being larger than the one-electron 
splittings caused by the environment, Hund's rules would still result in an m = triplet ground state, but with the 
lowest excited states also being spin-triplets.] 

Motion of the impurity away from the triangle center by a displacement (x, y) with ~x m. the direction of one of 
the three neighboring atoms, will give rise to mixing of the one-electron states via terms of the form 



n. 



split 



^1 X! \a\„aB,y{x + iy) + h.c. (15) 



-52 ^ \cL^R„asa{x + iy) + a\^asa{x - iy) + h.c 



<y=hl 



These will give rise to splittings of the doublet excited state that are linear in the displacement; this is just what- 
would-be the Jahn- Teller splitting if the doublet had been the ground state. If the displacement is directly towards or 
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Spin-triplet 
region 



Neighboring atom 



FIG. 1. Schematic representation of the minimum energy impurity positions (1,2,3) in the environment of the neighboring 
atoms (open circles). The solid triangle is the boundary between the regions in which the electronic ground state is a spin-triplet 
(interior) and a spin-singlet (exterior). Points 0, a, l3, 7 represent positions of the impurity at which the two lowest lying 
electronic spin-singlet levels are degenerate. At point these two states form an orbital doublet [see Fig.^ and Section 



[II 



1]. This figure corresponds to large ea, with the parameters defined in Section [II, 1 ea = 33. 3J = WK = 20(5i = 10(52 and 
units of X and y, such that a, jS and 7 are at distance approximately 0.6 from 0, see Fig.^. Depending on the choice of these 
parameters, the crossing points a, f3 and 7 may lie inside or outside of the triangle. Two types of lowest action tunnelling paths 
'A" and "B" are shown and discussed in Section III.l. 



away from one of the triangle corner atoms, say, with y = 0, the states can still be categorized as even or odd under 
reflection through that atom. The even state will also mix with the m = spin singlet excited state. 

We are interested in distortions large enough that the ground state is no longer the spin triplet, but one of the 
even or odd singlet states. It can readily be seen in simple models, like the one above, that this can occur, at least in 
principle. There will be then, for a classical impurity, three equivalent spin-singlet ground states, with displacements 
(i?, 0), {—R/2,±^/3R/2). In general, the lowest action tunnelling path between these positions will go around the 
excited state doublet degeneracy in the center of the triangle, and hence pick up a minus sign from Berry's phase. 

Unfortunately, there may also be other degeneracies in the spin singlet sector that can play a role. In particular, 
level repulsion between the even singlet states along the line y = (and two equivalent lines), can result in a second 
crossing between the even and odd singlet states. It should be stressed that these extra level crossings a,f3,j [see 
Fig.|l[ are distinct points, occuring only when the impurity lies on one of the three lines of symmetry, so that the 
two low-lying singlet states have different (even-odd) symmetry and hence can cross, see Fig.^j. In particular, if the 
modification of the effective interactions by the distortion is small, (e.g. for ea ^ J,K with no modification of the 
interactions that give rise to J and K due to the distortion), the main effect of the distortion will be the repulsion 
between the even singlet states; this situation is shown schematically in Fig||. In this case the degeneracies and region 
of spin singlet or triplet ground states are shown in Fig.^ If the lowest action tunnelling path is like that denoted by 
"A" in FigJ^, four degeneracies in the spin singlet sector will be encircled as the impurity moves around the triangle, 
resulting in four minus signs and hence a conventional symmetric orbital ground state of the tunnelling impurity, i.e. 
£ < 0. On the other hand, if the impurity tunnels via a path like that shown as "B" in Fig.|l|, only the triangle center 
degeneracy will be encircled and the tunnelling ground state will be the impurity L, R doublet. Either possibility 
can occur. It should be mentioned that if the repulsion between the electronic states (J, K) is also modified by the 
distortion, then the three extra degenerate points can be made to disappear in simple models, thereby leaving only the 
triangle center degeneracy and only paths of type "B" possible, hence making the impurity ground state a doublet. 

Note that spin flip effects can complicate matters somewhat, but will not change the overall conclusion. 

We have shown in this Section that the existence of an electronic degeneracy in either the ground state or some 
particular excited state when the impurity is located in a special, albeit energetically unfavorable, point in space (the 
triangle center) can be responsible for the restoration of that degeneracy for the combined impurity-plus-environment 
system. This only occurs in some regimes of atomic electron interactions and couplings to the environment. Clearly, a 
complicated calculation would be required to determine the amplitudes of the various processes reasonably accurately 
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FIG. 2. Energies of the low lying atomic electron states, discussed in Section IILl, as a function of the impurity displacement 
X along the j/ = axis [see FighJ. This plot corresponds to ea = 33. 3 J = IQK = 205i = 1052, the same as m FigJll When a; = 
the impurity is located at the triangle center and the electronic orbital doublet is between the '^5 spin-triplet, orbital-singlet 
and the spin- and orbital- singlet states. When x ^ Q (but y — 0)the orbital doublet is split linearly and all the states are 
now classified as even or odd under reflection through the x-axis. Due to repulsion between the spin-singlet even states, for 
a; < the lower even state crosses the odd one at point a as also shown in Figjl|. 



in an actual physical example. Furthermore, the above analysis is valid only when the impurity moves slowly compared 
to the atomic electron time scales. One should note that the above assumption—although intuitively obvious in 
many cases (especially away from electronic level crossings) is still controversial .E3 However the above discussion 
demonstrates that, at least in principle, the appearance of a Berry phase, and therefore an orbitally degenerate 
impurity ground state, is indeed possible. 



B. Effect of an excited impurity configuration 



In realistic situations, impurity tunnelling rates will probably be appreciable if the barriers between the three 
impurity positions are relatively small. Therefore the excited energies of intermediate impurity configurations that 
are comparable to these barriers, will be substantially lower than the bandwidth. As a result, a considerable fraction 
of conduction electron excitations will have energies large enough to allow transitions between the low lying and 
these excited impurity states. However, since the electronic excitation spectrum is continuous in the thermodynamic 
limit, the impurity motion cannot be treated adiabatically for these energies, as it was for the atomic electron states 
discussed above. Instead, the conduction electron assisted tunnelling of the impurity between the low lying and these 
excited states can be analyzed systematically by RG scaling. Here we show how the existence of an excited impurity 
configuration can lead to a positive effective E at low energies via coupling to conduction electrons. 

For simplicity, we only consider two sites, with the low lying impurity states at each site being d\. Let 
denote an intermediate excited impurity state, with extra energy U , that is large compared to the bare impurity 
tunnelling matrix element between sites 1 and 2, but small compared to the bandwidth W {U/W <C 1), see FigJs. 
In this regime we can perform a two-stage renormalization procedure: Initially we include the excited state, allowing 
electron assisted tunnelling events between states 1, 2 and E, these may have tunnelling amplitudes ti much greater 
than £. Then, when the renormalized bandwidth becomes of the order of U or smaller, excitations of the impurity to 
state E will become negligible, in which case we can drop the excited state from the effective Hamiltonian. Thus the 
effect of the impurity excited state can be described by the following terms 

An= Ud^dE +yJ2 {AaC2a + 4.^1,) 

(7 

+2TTti ^ f d^dic|^c_E<x + d^^d2c\^CEa + h.c.\ (16) 
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where c^^ is an extra local electronic state, which is even under the exchange of site indexes 1 and 2 and y is the 
average mixing amplitude of conduction electron states 1 and 2 in the presence of the impurity. Note that in the 
absence of the impurity, y (together with the other t^ms in A7i) will vanish, restoring the orthogonality of states 1 
and 2. In general, both ti and y are energy dependent .E3 For simplicity, we only consider their constant part, evaluated 
at the Fermi energy. Ener-gy dependent deviations will introduce corrections but will not qualitatively change the 
behavior describpi here.E3 A crucial observation is that y can be positive or negative due to two effects. First, as it 
was seen in Ref.Ea, since c\^, involve electronic excitations in the conduction band, y is an oscillatory function 
of kpR, where kp and R are the Fermi wavevector and the intersite distance, respectively. Roughly, this oscillatory 
behavior corresponds to the amplitude of an electron at the Fermi energy tunnelling between sites 1 and 2. Second, 
since it is the presence of the impurity that causes the mixing of the otherwise orthogonal states 1 and 2, y shouM 
also depend on the impurity-electron coupling. As a result, its sign should also depend on the sign of the coupling.tS 




FIG. 3. Schematic illustration of a double well potential for an impurity with low lying states d\, d\ (and their representative 
wavepackets localized around sites 1 and 2, respectively). The wavefunction of the excited impurity state at an energy U 
higher than d\, d\ [see Section 



[II 



2], as shown, has significant amplitude in both wells. 



Neglecting all other terms involving the impurity, in this regime the RG equation of £ normalized by the conduction 
bandwidth W is 



— =£ 
dl 



6yt' 



(17) 



while y and ti are marginal (for weak electron- impurity interaction). The second term ceases to contribute to 
the renormalization of £ once the renormahzed bandwidth becomes of the order of U. At that scale e' — W/U, 
£{l) ~ {£ -\-Qyt\)W /U so that (after reinstating the renormahzed bandwidth factors) the effective renormahzed £ that 
should appear in the lower energy effective Hamiltonian, Eq(^, is 



£f=i£ + GytjW, 



(18) 



independent of U. Below the energy scale U, the impurity is limited to hops between the low lying states d\, d\ (and 
d\), with an effective hopping £. Note that since the state has appreciable amplitude in both wells, t\W may be 
larger than £. Therefore yt\W may be comparable to £, and, since the sign of y is arbitrary, £ may be positive or 
negative, whatever £ may have been originally. [Note that there will also be contributions to £ from hops 1 ^ E ^ 2 
not involving electrons, but these need not dominate over the electron-assisted process included here. The effects of 
an excited impurity configuration were first considered by Zarand and ZawadowskiE3 in the context of an impurity 
tunnelling between two sites, in order to show that the one-electron assisted tunnelling between the low- lying levels 
(which, for simplicity, we did not consider here) can be enhanced. However, they did not consider the equally strong 
enhancement of the tunnelling rate of the impurity alone, which we analyzed above.] 

In summary, we have presented three situations which can result in an impurity L,R doublet ground state at low 
energies. The first two (Section pl|.l) involved the effect of localized atomic electrons with energies far from the Fermi 
level. In the third case (Section p^H|.l) we saw that, closer to the Fermi level, conduction electrons in the presence 
of excited impurity configurations can renormalize the effective bare tunnelling element £, so that it can take either 
sign. After all these effects have been taken into account, at intermediate energy scales the coupling of the remaining 
conduction electron degrees of freedom to the impurity can be represented by the effective renormahzed coupling 
terms introduced in Section |l| these will now be analyzed. 



9 



IV. WEAK TUNNELLING CROSSOVER 



In most physical situations, it is likely that the bare impurity tunnelling £ and the electron assisted tunnelling 
processes will be small. If the impurity-position-dependent electron-impurity interactions Vi and V2 are also small, 
then there will be a simple crossover at an energy scale of order £. At lower energy scales, the processes that involve 
ds will be strongly suppressed; those remaining will be Ai, Qi and A3. By this crossover scale, the various Qj and 
Aj, as they are all marginal for weak coupling, will only have been modified from their original values by amounts 
quadratic in the various couplings; £ will have been renormalized similarly in addition to its simple growth under 
renormalization with eigenvalue one. For weak coupling, two-electron terms, like those in nf^^p, will be irrelevant with 
eigenvalue —1 and hence will be negligible at the crossover scale. 

Beyond this crossover scale for weak interactions, the couplings other than Ai, A3 and Qi can be ignored (at least 
for now), but these will renormalize each other. With all three of these small the flow equations in this regime are 

dQl nA2 



dl 
dAi 



dl 

(reflecting the approximate SU{2) symmetry for Ai = Qi) and 



2Af (19) 
2Q1A1 



dl 



Q1A3. (20) 



If Qi < — I All, then the system is like a ferromagnetic Kondo problem so that Ai flows to zero and the ground 
state will be an impurity orbital doublet. We thus restrict consideration to Qi > — |Ai|. Since Qi involves impurity 
hopping but Ai does not, we expect \Qi \ < |Ai|, although this need not be the case. 

We thus see that A3(/) will tend to decrease in the interesting regime, while Qi and Ai increase. The Kondo 
temperature Tk, below which the intermediate coupling behavior will obtain, is determined by the energy scale at 
which Ai(Z) becomes of order unity, yielding 



Tk exp 



|Ai| 



(21) 



with W the electron bandwidth and k depending on the ratio of Qi and Ai initially, or more precisely, at the crossover 
scale at which £{l) is of the order the renormalized bandwidth; specifically for \Qi \ < |Ai|, 



1 

^=^= arctan , 



2 Wl 

% (22) 



V 

which is 7r/4 in the small Qi limit. 

If the interactions of electrons with the impurity are reasonably strong, i.e Vi and/or V2 are not small, then the RG 
flow away from weak tunnelling is more complicated. In particular, the weak tunnelling eigenvalue of the impurity 
hopping, Ag is modified to 

X, = 1-^{V,' + 2V,') (23) 

by orthogonality catastrophe effects. This expression is valid for Vi and V2 of order unity if they are replaced by the 
associated phase shift-like quantities, see Ref.113. The renormalization of the various electron-assisted tunnelling terms 
Ai, . . .Aq will also be modified by the impurity-electron interactions, as discussed in Ref.E£l. Indeed, for sufficiently 
large Vi and/or V2 some of these can be more relevant than £, although in this case the feedback of these into £ will 
become important. Likewise, £ combined with Timix will generate the electron- assisted hopping terms if these do not 
already exist. 

Thus for strong interactions, the weak tunnelling crossover will be rather complicated with several processes poten- 
tially playing significant roles. But the coupling of primary importance for the Kondo-like physics, Ai, can itself be 
quite large initially; in this case the Kondo temperature will essentially be when £ becomes of order the renormalized 
bandwidth i.e. 
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(24) 



with Ag < 1. Note that initially Qi will not be generated by the Vi terms, although they will affect its renormalization, 
but once the effective £ becomes appreciable, Qi will be generated by Ai as in Eq(19) and then rapidly become of 
order unity. 



A. Symmetry Breaking Processes 

The dominant C^^ symmetry breaking process is given by Eq(p^. For weak coupling M is relevant with eigenvalue 
one. Beyond the crossover energy scale at which £{l) becomes large, M can be renormalized non-trivially by Qi and 
Ai. If Ai is large, there will not be a substantial regime in which this will occur and we can simply use the magnitude 
of M at this crossover as an input to the flows in the interesting intermediate coupling regime analyzed in the next 
Section. On the other hand if Ai (and Qi) is small, then there will be a regime of non-trivial renormalization of M, 
with corrections to its eigenvalue of order A^ and Ql . But since the logarithmic range of scales over which this will 
obtain is of order I ^ 1/Ai, the net effect of this will be small. 

Thus in both regimes, the weak tunnelling crossover will not appreciably alter the effective magnitude of the 
dominant symmetry breaking perturbations. 



V. INTERMEDIATE COUPLING 2CK FIXED POINT 



We now analyze the intermediate coupling regime and show how two-channel Kondo-like behavior arises. We work 
in the regime where the gap between the symmetric impurity state and the lower lying Fa-doublet is much larger than 
the renormalized electronic bandwidth, i.e. at energy scales much less than the Kondo temperature Tk- In this case 
the impurity is essentially frozen in its doublet ground state and we can therefore impose the constraint 

dldL + d^j^dn = 1 (25) 

and limit the analysis to processes involving only these two states. In this Section we first identify the 2CK fixed 
point, which is located at intermediate values of the important coupling parameters Qi and Ai and therefore cannot 
be accessed using a weak coupling analysis. We then show that it is stable when C^y is an exact symmetry. Finally 
we find the set of leading irrelevant operators which will control the low temperature behavior of the system. 
It is convenient to introduce pseudo-spin operators for the impurity doublet 

a, = dldL - d^iidR (26) 

and 

a+ = dldn (27) 
cr_ — d^j^dL 

In order to proceed with our analysis we bosonize the electronic degrees of freedom by introducing bosonic fields 
(x) for fi = L, R, S and a , J, in the following way: 





(x) 


1 




CR] 


(x) 


1 




CST 


(x) 


1 




CLi 


(x) 


1 




CRl 


(x) 


1 

y/2'KTc 


g'*J?4(a;)gi7r(JVT+Afi) 


csi 


(x) 


1 

\J2-KTc 


g«*Si(a:)gi7rJVsigi7rJV 
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where the exponentials with the number operators 



N, 



{.La 



dx C' {x) Cf,a{x) 



(29) 



have been inserted to insure anticommutation relations between the various different fermion operators. Also, 

is an cutoff energy scale of the order the renormalized bandwidth. In the above notation the (normal ordered) local 

electron operators can be expressed simply in terms of the bosonic variables 



'^liaix) Cfj_a-{x) — 



J_d%Ax} 
2tt dx 



(30) 



The electronic kinetic energy can also be represented in a standard manner in terms of bosonic variables (e) which 
obey the canonical commutation relations 



with 



by 



^^Lcr (x) 



de 

/27re 



f°° de 



(31) 



(32) 



(33) 



with both integrals involving only positive energies (corresponding to the fact that the electrons cj^g. (e) are only left 
moving) and with a cutoff energy of the order of the renormalized bandwidth t~^. [Note that the energy ^-function 
in Eq(^l|) is rounded at energies e « t~^.] It is useful to define even and odd combinations of the L, R Bose fields 



V2 



Then the most important parts of the effective Hamiltonian can be written as 



Atn H T=C' 



^ ^ dx 



(34) 



(35) 



Ai 



E 



{aj^i]a exp [^1^2$ OCT + h.c^j 



where the boson fields in the last two terms are evaluated at the origin (x — 0) and we consider later the term 
proportional to A3, which will turn out to be irrelevant. Here we have defined the factors 



'la c , 



(36) 



which can take values ±1, since the number operators A^^o- can only take integer values. In the case that t]^, has a 
definite sign, (i.e. when TV^i + 7V/{| and Nj,^ + Nb.i are separately conserved) these factors can be neglected, as we 
shall do for the moment. Note that at this point the S'-electrons are decoupled; we will later reintroduce terms that 
couple them to the L and R degrees of freedom. 

Eq(p3) is essentially equivalent to the Hamiltonian of the anisotropic 2CK model, which has an 0(2) symmetry 



^/V2, 



f7± 



SU(2) symmetric 2CK fixed point.QMtZI 

After decomposing the bosonic fields and r\a into charge and spin components. 



-fTx]; this is sufficient to make it flow to the 



cr*es) 



etc. 



(37) 
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and 



'/c — c I 



(38) 



with a — ± for spin |, |, respectively, we follow Emery and Kivelson,l3 and perform a unitary transformation 



U = exp 



which transforms Hef / to the for 



H'eff — UHeffU^ =1-10 + 



2Ai 



(Tx COS$os 



(39) 



(40) 



As was first pointed out by Emery and KivelsonO the above Hamiltonian can be diagonalized exactly when — 1/2 
via refermionizing the bosonic and pseudo-spin operators. This solution is analogous to the ToulouseJimitEa of the 
one-channel Kondo model; here it captures the behavior of the 2CK model. Sengupta and GeorgesEil showed that 
az ^g"" is ona of the leading irrelevant operators with dimension 3/2, while a second dimension 3/2 operator was 
found by usEZi in the context of the 2CK model, consistent with expectations from conformal field theoryjjd The 
RG eigenvalues here are one minus the dimension of the local operators. There are no relevant operators consistent 
with 0{2) [and spin SU(2)] symmetry. However, the underlying symmetry of our system is C'^y, which is lower than 
0{2). Although no extra relevant operators exist if we lower the symmetry to C^y, there will turn out to be one extra 
irrelevant, dimension 3/2 operator allowed and a fourth, if electrons from the antisymmetric representation r2 are 
coupled to the impurity. 

Unfortunately, the refermionization technique does not work for all important operators allowed by C^y. Thus, in 
order to show that the above are indeed the onbi. allowed leading irrelevant operators, we analyze the problem with 
a method recently T[Lsed by FabriziQ_and GogolinEj in the context of the four-channel Kondo model following ideas of 
Fisher and ZwergertZl and Schmid.tJ 



A. Analysis of effective Hamiltonian 

The approach described below is based on the observation that Ai in Eq(pO[) formally has dimension 1/2 and 
therefore will flow under RG to arbitrarily large values. In this case it makes sense to "diagonalize" this term first 
and treat Tig as a perturbation. Without loss of generality we assume that Ai > 0. Neglecting TYq, the minimum 
energy values of <I>os {x — 0) are, in the impurity basis where ax is diagonal, ^os {x — 0) = (2n + 1) vr when = +1 



and <I>os (a; = 0) = 2n7r when ax — —1- The presence of Ho has two effects. First, it includes (via Eqs 32 and |33|) 



the conjugate momentum of <I>os (0), Hos (0) = Thus, in the semi-classical approximation, which is valid 

for large Ai, the phase (0) will tunnel occasionally between its various minimum energy values, with an effective 
tunnelling-rate proportional to e"'^*'"'' , where Skink is the classical action of the motion of $os (0) in the inverted 
potcntialo (Ai — Ai), from, say — tt to tt. The value of Skink, which is large_svhen Ai is large, is not important 
and clearly non-universal, since it depends on the details near the cutoff energy.c3 

The other effect of Tin is via the coupling of {x — 0) to bosons located away from the impurity [^os {x ^ 0)], 
which will tend to screen the "kinks" (tunnelling events) of <i>os (0). The underlying reason is that the ground 
states of $os {x ^ 0) with different values of $os (0) are mutually orthogonal, the manifestation of Anderson's electron 
orthogonality catastrophe. To analyze this effect, we integrate out the bosons $os (x 7^ 0), leading to the following 
effective imaginary time action for [x = 0, r) = $os (t) 

{27:f J J {r-T'f ^ ' 

with $os (t) typically being a series of widely separated (on the cutoff scale) kinks betwepi different minima. This 
action has been studied extensively in a variety of quantum impurity problem contextsli3E3c3C3 and it is straight- 
forward to show that the operators Cm, corresponding to $os (''') hopping by m7r, have dimension Therefore, 
when ax has no dynamics, the leading irrelevant operator is 02, corresponding to hops between adjacent minima of 
the cos$„, term. 
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B. Leading and other Irrelevant Operators 



In the absence of relevant operators, at temperatures small compared to Tk the system will flow towards the fixed 
point. The flow close to the fixed point will be dominated by the leading irrelevant operators, i.e. the ones with the 
least negative RG eigenvalue, since other operators, being more irrelevant, will renormalize to zero faster and therefore 
will become negligible. At non-zero temperatures the RG flows will be terminated when the renormalized bandwidth 
becomes of order the temperature. At this scale the only non-negligible processes that will enter in physical quantities 
are the leading irrelevant ones and the temperature dependence of their renormalized amplitudes will determine the 
low temperature behavior of the system. 

First we analyze the leading irrelevant operators with dimension 3/2. As mentioned above, the first is azd^ocl dx. 
Since <^oc is decoupled from and ax when Qi — 1/2, the derivative term clearly has dimension 1. The effect of ct^ 
(orj^quivalently dy) is non-trivial: since it anticommutes with Ux in the Ai term, it effectively reverses the sign of 
tTa; Hj As a result, the ground state value of shifts by tt when is applied, thereby giving (and ay) dimension 
1/2, since it forces a tt kink in <i>os, just like Oi above_As a result, the operator proportional to (Qi — 1/2) has overall 
dimension 3/2, in agreieHient with other approaches oIlJ 

As discussed in Ref.EZI, there is a second dimension 3/2 operator in the spin sector, which in the original electron 
operator representation is 

AdR 4|CLx |c[^CLT + 44CRt I + c|j^CL| {c^^cli + c\^cri | (42) 



h.c. 



written in an explicitly spin-S'C/ (2) invariant form. Note that this operator involves two electron processes [part of 
^iop] ^^'^ thus is irrelevant for weak coupling. After bosonizing and performing the unitary transformation of Eq(^9|) , 
the most relevant part in the square brackets becomes 

• ^ ^^^^^ lA-i^ 

o-ysm^os-^— . (43) 

Apart from the above operators there is an extra leading irrelevant operator consistent with the C^^ symmetry. 
Similar to the above, in the original fcrmion representation, this operator has the spin-S'C/ (2) invariant form 



+ h.c. . (44) 



After bosonizing and performing the unitary transformation of Eq(p9|), the last term in the square brackets becomes 
proportional to 

gySin^os ■ (45) 

exactly analogous to Eq(^). Since ay flips the sign of ax in the Ai-term and is essentially frozen, Eqs(^5|) and 
( ^ ) have dimension 3/2. Note that although nominally sin$os = Q at the minima of the cosine potential, the kink 
forced by the ay gives sin $os a non-zero value at the kink. 

It is important to comment that the extra dimension 3/2 operator (Eq(p5|)) only exists because of the presence of 
the symmetric electron state c^^., which is invariant under C^y. In fact, any C^v singlet fermion bilinear made up 
of electrons that can couple to the impurity through its spin, like cs here, will produce a dimension 3/2 operator. 
In particular the operator similar to Eq(^^ , but with the cg-electrons replaced by the ca electrons of the odd one- 
dimensional r2-representation, will also have dimension 3/2. Note however, that due to their symmetry, r2-electrons 
tend to interact more weakly with the impurity. 

The first two terms in Eq(^), transformed via Eq(|39|), become proportional to 

ax cos <^es - V2<i>ss (46) 

where we have neglected the exponential factors exp [i {N^c ± Nes)]- Naively this operator also has dimension 3/2. 
However, due to the presence of $es, we need to be more careful with the exp [i {Nec ± Nes)] in the Ai term of Ti-'eff^ 
Eq(^). Including these factors explicitly alters the dominant terms in Ti'eff to 
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2A 

Hef / = Ho H Wx COS TrA^ec + Cfy sin -nNec) cos [$os - 7I'A'"es] (47) 

Thus the presence of an exp ±i$es factor in Eq(p6[) causes a sign change of the Ai term since it takes N^s — > N^.^ + 1 
and therefore behaves just as CTz does, increasing the dimension of the operator to 2. Similar arguments hold for 
subleading terms from Eq(^2|). 

Naively one would expect all parts of operators such as Eqs(^) and ( ^^ to contribute to the leading irrelevant 
operators since the they are orbital and spin singlets. The reason they do not is that through the bosonization scheme 
implemented above, in addition to the unitary transformation of Eq(p3), the underlying orbital SU{2) symmetry of the 
Hamiltonian is broken down to an 0(2) symmetryllZI (which is the symmetry of the XXZ anisotropic orbital Kondo 
modelO). Therefore some parts of the various operators will mix will orbital-spin-2 operators associated with this 
symmetry breaking, thereby acquiring different scaling dimensions. Since the orbital-spin-2 operators have dimension 
2 at the orbitally SU{2) invariant fixed point and are thus more irrelevant than the dimension 3/2 leading irrelevant 
operators, the breakdown of the orbital SU{2) symmetry to an 0(2) symmetry will not affect the low temperature 
behavior, as the dimension 2 operators will become negligibly small at low temperatures. This is the underlying 
reason why the solvable point of the 2CK mndol found by Emery and KivelsoiC captures correctly the singular low 
temperature behavior of the 2CK fixed point .ljEJ 

Other than Ai and Qi, the only one-electron impurity operator that survives when ds is suppressed, is A3. The 
most relevant part of this near the 2CK fixed point can be written as 



V2 2 



V2 2 



(48) 



Naively this has dimension 15/8, since $os is frozen. However, the presence of expii^^ shifts Nes in Eq(^7|) by 
±1/2, turning the cosine into a sine. This forces a shift of by ±7r/2, thereby increasing the dimension of the A3 
operator in Eq(^) by 1/8 to 2. This is in complete agreement with the dimension of A3 derived using conformal 
field theory methods. To see this, we observe that this operator combines csa or c^^, which are free Fermi fields, 
uncoupled from the impurity at the fixed point, and hence having dimension 1/2, to {^\dj^c\^ + d^j^d^c^j^^ and 

{^\d]iCRrj +d^j^dLCL^, respectively. These operators by themselves would break the C/(l)-charge symmetry of the 

model; since they involve an odd number of electron operators, they have charge 1. [In our bosonization scheme, this is 
manifested in them changing sign under <I>ec $ec + 27r]. Furthermore, they clearly transform as a spin-doublet. The 
dimension of a charge-1, spin-1/2 operator at the 2CK fixed point is n -f 1/2, with n a non-negative integer .123 However 
the operators with n = do not respect the C^ri symmetry, in particular, unlike the above operators, they are not 
invariant under cl exp and cb, cr exp [— ^] , since they transform as doublets under 0*3^. As a result, 

the allowed operators have n = 1 so that the scaling dimension of {^\dRc\^ -f d^j^dLC^j^^ and {^d^i^djiCRa -\- d^j^dLCLaj 

is 3/2; combined with the 1/2 of csa or Cg^, this makes the dimension of A3 equal to 2, as found above. 

By similar arguments, it can be shown that the dimensions of all physical operators are integers or half integers. 
This should be expected from the conformal field theory point of view, since the 2CK fixed point only has boundary 
operators with half-integer dimension.E3 Combining these with free fermions, such as csa, that couple only through 
their spin (which plays the role of the "flavor" in the spin 2CK model), does not alter this result, as we saw explicitly 
in the case of A3. 



C. Physical consequences 



The low temperature behavior of the impurity and the electrons that couple to it will be governed by the 2CK fixed 
point discussed above and the flow towards the fixed point. In the absence of the irrelevant operators, the effective 
Hamiltonian of Eq(|4^) with Qi = 1/2, i.e. Eq(p7|), couples only the ^tj—Bose fields to the impurity. [Although strictly 
speaking this is not the fixed point Hamiltonian, as discussed in RefJlJ it has the same essential structure, differing 
from the fixed point primarily in Ai being finite.] 

This means that quantities like the spin susceptibility, x, which involve only other linear combinations of the Bose 
fields-in this case correlations of the spin operators from electrons in different representations, e.g. ^^g^*-"-* -will 
decouple from the impurity and be non-singular. However the leading irrelevant operator in the spin sector Eq(|4^) 
can combine with to yield CTj, sin <I>os , which has dimension 1/2 as discussed above. Similarly, the other two 

leading irrelwant operators involving ^'-electrons and A-electrons can also yield (T^ sin $os. As shown by Sengupta 
and Georges ,Eil this effect gives rise to a contribution to the susceptibility 
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Ax ~ In (Tk/T) 



(49) 



at low temperatures, with the coefficient being proportional to the sum of the squares of the coefficients of these three 
leading irrelevant operators in the effective low energy Hamiltonian. 

On the other hand, all dimension 3/2 operatopSi, including the one proportional to (Qi — 1/2), will contribute to 
the specific heat Cy, giving a singular correctionE^ 

ACv -Tlii{TK/T) (50) 

with the coefficient being proportional to the sum of squares of all four leading irrelevant operators. As a result, 
the normalized ratio of TAx to ACy-the Wilson ratio-will not be universal. However, due to the similarity in their 
form, the three operators entering the susceptibility, §f fj, sin$os, with the $ in the derivative being one of <f>es, $5s 
and $As, will enter the specific heat in exactly the same combination. Thus for Qi ~ 1/2 the normalized ratio of 
the specific heat over the susceptibility will be universal. As a consequence, to find a general universal ratio, one 
extra independent measurable quantity is recwked, that involves, e.g. only (Qi — 1/2)^. It should be noted that the 
contribution to the resistivity which scales ag3'Q 

Ap^Vf (51) 
is not suitable, because it involves a sum over the amplitudes of all the leading irrelevant operators. 



D. Symmetry breaking perturbations 

The most important symmetry breaking perturbation is a strain or other anisotropy that favors one of the three 
impurity positions. This most simply appears in the Hamiltonian as in Eq(p^) where the phase of M determines the 
preferred direction for the impurity relative to the center of symmetry. This corresponds to axReM + ayImM in 
the effective Hamiltonian, operators which have dimension 1/2 and are hence relevant, albeit less so than from naive 



power counting. At energy scales above Tk, M just grows linearly with energy [as described in Section IV.l], but on 
scales below Tr-, it grows more slowly with eigenvalue 1/2. The perturbation thus becomes important and destroys 
the non-Fermi-liquid behavior below the energy scale 

T„~^, (52, 

the |M|^ reflecting the eigenvalue of 1/2. Thus the novel 2CK behavior is more stable than might be expected. Note 
that similar crossovers are found in the presence of spin and electron- orbitally-coupled magnetic fields since these 
break the spin symmetry or the L, R symmetry, respectively. 

It should be noted that the operator M in Eq(p^) can also be generated through electronic processes that break the 
triangular symmetry caused, e.g., by the presence of static impurities in the vicinity of the tunnelling impurity. These 
give rise to mixing between the electronic wavefunctions interacting with the tunnelling impurity or between these 
and other wavefunctions.Lj This mixing breaks the Csy symmetry. In the presence of electron-assisted-tunnelling 
processes, such as ^i, . . . , discussed before Eq(P), these C^y symmetry breaking electronic terms, such as 



E 



PcLcna+P*cl^cL^ (53) 



will generate under RG the operator M. 



E. Effects of higher symmetries 



Impurities with a Fs-doublet ground state tunnelling around a point in the lattice with higher symmetry, such as 
octahedral symmetry in a cubic lattice, will clearly not have any relevant operators at the 2CK fixed point. The 
reason is that since these relevant operators are not allowed by C^y, a subgroup of the octahedral group O, they 
will not be allowed by O either. However, there will be extra irrelevant operators, which turn out also to have 
dimension 3/2, therefore not essentially changing the low temperature physics of the 2CK fixed point. Electrons 
from representations T such that F (g) F contains the trivial representation, will also have an associated dimension 3/2 
operator of the form of Eq(^). Furthermore, since F4 (g) F4 and F5 (g) F5 contain, apart from Fi, also the doublet F3, 
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electrons in these representations can couple orbitally with the impurity, in the form dj^dji Cp.^crio-, for i = 4,5, 
giving two additional dimension 3/2 operators. As a result, in the case of octahedral symmetry O, where there are 
two 3-dimensional representations of electrons, r4 and Fs, there will be four extra dimension 3/2 operators, leading 
to a total of eight leading irrelevant operators. 

If the impurity is in a triplet ground state or r5 in the octahedral, or T in the tctrahedral case), due to arguments 
similar to those in this paper, it is anticipated that there will be a stable non-Fcrmi-liquid fixed point, analogous to, 
but different from, the 2CK fixed point for an impurity Fa-doublet ground state. In this case the situation is similar 
to thc_critical point of the doublet case at which £ changes sign, corresponding roughly to a triplet impurity ground 
state;o this is briefiy discussed in the next Section. 

VI. CRITICAL MANIFOLD SEPARATING THE 2CK FROM THE FERMI LIQUID FIXED POINTS 

For all of the physics discussed so far, the sign of £ is crucial. Here we consider what happens if £ can be made to 
change sign by a change in the system, for example by applying pressure. From the analysis of Section it should 
be clear that it is really the sign of the effective hopping, £, which includes renormalizations due to the effects of the 
electron assisted tunnelling, etc., that is important. Generally there will be a critical manifold separating the novel 
2CK low temperature behavior from a conventional Fermi liquid regime in which the impurity is in its non-degenerate 
symmetric state with simply potential scattering of the electrons. This is schematically shown in Fig.^. Since both the 
fixed points are stable at zero temperature, the critical manifold separating their domains of attraction can generically 
be crossed by a one parameter change in the system. 

When the system lies on the critical manifold, which corresponds to £ = at small coupling, it will presumably flow 
to a new unstable fixed point. We conjecture that, apart from the charge-?7(l) and spin-5'C/(2) symmetries, this has 
an additional orbital-5'C/(3) symmetry, where the L, R, S triplet is identified with the fundamental representation pf 
SU{3). The behavior of the system on the critical manifold and close to this fixed point will be analyzed elsewhere,E2l 
but below we quote some of the important results. 

Within the critical manifold this fixed point is stable and the RG fiows close to the fixed point are determined 
from the dimension of the leading irrelevant operator, which in this case is 7/5. Therefore the impurity corrections 
to the specific heat, resistivity and susceptibility at low temperatures are ACv T^/s, Ap ~ r2/5 and Ax ~ T-i/^, 
respectively. Note that although the leading-jirrelevant operator and the symmetry of this fixed point are the same 
as those of the three-channel Kondo model,cil the fixed point itself and the operator content are not the same. This 
can easily be seen by noting that while in the three-channel Kondo model an impurity S'[/(2)-doublet couples to 
3 channels of electrons through their spin, in this case the impurity is an SU (3)-triplet coupling to 2 channels of 
electrons through their orbital SU (3) index. 

Deviations away from the critical manifold are relevant dimension 3/5. In a C^y symmetric situation only one 
such operator exists, which can be represented by {2dgds — d\dL — d^d_R)/3. This is therefore analogous to £ and 
it governs the RG flows away from the critical fixed point toward the Fermi-liquid and the 2CK stable fixed points. 
Therefore, tuning only one parameter (such as external pressure) the novel non-Fermi-liquid fixed point associated 
with this higher orbital symmetry might be observed. In particular, close to the critical pressure, there will be a low 
temperature regime with the critical behavior crossing over below a temperature 

Tx^\p-Pcf^ (54) 

to one of the two possible low temperature behaviors. 

It should be noted that the above behavior also holds for the case of spinless impurities with a triplet degenerate 
ground state in tetrahedral or octahedral sites in a cubic crystal. In that case there will be no unstable directions 
close to the fixed point in the absence of symmetry breaking terms, just as for the 2CK fixed point in the case of Cst,, 
and the low temperature physics would be governed by the critical fixed point and its corresponding non-Fermi-liquid 
behavior. 



VII. CONCLUSIONS AND RELEVANCE TO EXPERIMENTS 

We finally discuss briefly the possible relevance of impurity tunnelling to the nanoconstriction experiments of Ralph 
and Buhoa an ,113^113 as well as potential future experiments. 

In Rcf.EZi we showed that the inteijpretation of the observed non-Fermi-liquid-like behavior of these experiments as 
a defect tunnelling between two siteaU runs into trouble due to the existence of two relevant operators, that drive the 
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FIG. 4. Schematic of renormalization group flows for the three site problem at zero temperature. There are two stable fixed 
points, namely the non-Fermi-liquid 2CK fixed point and the conventional Fermi-liquid fixed point, at which the impurity 
ground state is a Fi-singlet which is completely screened at low temperatures. The critical manifold separating the basins of 
attraction of the two stable fixed points is also shown. The flow on this manifold goes to a critical fixed point, which has higher 
orbital symmetry. The critical fixed point has one unstable direction that respects Ca^; the dashed lines show the universal 
crossover fiows in this unstable direction from the unstable to the stable fixed points. 



system away from the 2CK fixed point and therefore make it highly improbable that a typical impurity could lie close 
enough to the critical manifold for the non-Fermi-liquid behavior to be observable. Another possible interpretation 
is in terms of an interstitial, vacancy or other impurity in the copper lattice structure which has a Fa-degeneracy. In 
such a situation, the relevant perturbations would only be induced by strains or other such effects; these tend to split 
the doublet, as discussed above. r~\\~\ 

In the copper nanoconstriction experimental '^rl the apparent Kondo temperatures are quite large, approximately 
Tk ~ 10— 15-/^. This suggests that the impurity-electron interactions are appreciable and thus the Kondo temperature, 
within our model, would be the energy scale at which £ becomes large. [As seen in Section the RG eigenvalue of 
8 should be less than one. However, for not too large Vi and V2 a lower bound for £ is Tk-] Such a fast tunnelling 
rate suggests that a hopping process might involve either the slight rearrangement of several atoms, as this can have a 
lower effective mass, or a small displacemenLof an impurity (or impurity complex) from its central symmetric position. 

In the copper nanoconstriction devicesEil no significant deviation from the y/T temperature dependence of the 



impurity conductance signal has been seen down to 50mK, the lowest temperature studied. From Eq(52) one can 
therefore estimate that the splitting epergy scale \M\ would have to be less than 0.7K, i.e. less than 0.05£. In 
more recent experiments in Pd and VlB a larger crossover temperature of about 1.5K was observed, which would 
correspond to M ~ 5K, which is about O.STr-. The stronger splitting in these materials-piight be attributed to the 
fact that the nanoconstrictions are more stressed than the corresponding copper devicesla. 

In order to see whether a scenario involving impurities with triangular symmetry might be reasonable, we need a 
rough estimate of the magnitude of the effects of strain. To be concrete, we assume a Jahn- Teller potential, which in 
polar coordinates about the center of the triangle has the form 

Vjrir) = -r^ - Br^ cos 361 + Cr^ (55) 
- Ar - AV^ cos 36* - e • £> • r 

with e the strain tensor and A and A' respectively the linear and quadratic splittings due to the Jahn- Teller effect. 
[Note that the potential for a Jahn- Teller impurity with cubic symmetry that has a doubly degenerate ground state 
would have a similar form to Eq(|55|) . In this case r would be a tensor quantity, even under inversion, that parameterizes 
the deviations of atoms close to the impurity from their equilibrium positions.] If r is measured in units of the lattice 
constant, a, then A, B, C and the third rank tensor D are all of order elastic energies IQ^K. In the absence of strain, 
there will be three minimum energy positions of the impurity with, for A and A' ^ A/a small, r « rg ~ A/A. These 
will be separated by barriers of order (ro/a) 10^ K, with comparable contributions arising from both the second and 
fifth terms in Eq(p5|). The tunnelling rate will be 
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£^ujD\—e'^ (56) 
V a 

with LUjy the Debye frequency and S the tunneUing action. Thus for £ to be about lOK, the action S cannot be very 
large. 

The dominant effects of strain arise in two ways: first, the tunnelling rates between the three sites can become 
asymmetric. Since some of the barriers will be changed by of order eDro from Eq(|55|), we see that the effective M 
resulting from this will be 

MtunnelHng ^ 5£ ^ ° ^ £e [ — ] (57) 

which decreases with rg. In contrast, the direct splitting from the change in energy of the minima is 

MspHt ~ eDro - e (^) lO^K, (58) 

which grows with tq. For both of these to be less than IK, one would need rg ~ O.lA and e ^ 10~^. So small a strain 
appears implausible in a nanoconstriction. Thus, unless an impurity configuration with anomalously weak..cQupling 
to strain or with other special properties were to exist, this explanation for the experiments of Ralph et atJtJ seems 
unlikely. 

Otlier possible interpretations exist, in particular, a possible explanation of the data was proposed by Wingreen 
et a£3 that involves the existence of static randomness in the neighborhood of the constriction. But overall, the 
interpretation of these experiments is controversial at this time. 

But impurity tunnelling about—symmetric positions in metals might occur in other situations. For example, it 
has been known for some timeEju that undersized impurities often form rather stable complexes with interstitials. 
These complexes, in fee lattices, form octahedral configurations. Thermally activated hopping between equivalent 
configurations ofr-these complexes has been experimentally demonstrated at low temperatures in the case of dilute 
CoxAli-x alloys O^hich suggests that low barriers might exist, even though in this system the tunnelling rates are 
probably very lowB3 

In conclusion, we have shown that an impurity in a metal tunnelling between positions with high symmetry~in 
particular, triangular symmetry-can lead to generic non-Fermi-liquid behavior, which can be characterized by the 
two-channel Kondo model. Symmetry breaking perturbations were shown to be partially screened, thereby reducing 
their effects and widening the potential temperature range for observation of the novel 2CK behavior. In addition, 
we have shown that by tuning a single parameter, such as pressure, the impurity behavior can change to that of 
a conventional Fermi liquid, but at the critical pressure the behavior would be controlled by a fixed point that is 
conjectured to have a higher SU (3) orbital symmeitijy. A detailed analysis of the behavior on and close to this critical 
manifold will be discussed in a future publication.E3 
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